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Preface  and  Abstract 


The  main  important  difference  between  the  dynamics  of  the 
atmosphere  and  of  the  oceans  is  concerned  with  the  driving  forces. 

The  atmospheric  circulation  may  essentially  be  compared  with  a 
thermodynamic  engine  driven  by  the  temperature  differences  between 
low  and  high  latitudes  on  the  rotating  earth,  and  this  circulation 
is  to  some  degree  modified  by  the  land  and  water  distribution  and 
by  oceanographical  conditions  at  the  sea  surface. 

In  the  oceans,  the  temperature  (and  salinity)  differences  be- 
tween the  equator  and  the  poles  play  only  a minor  role  as  long  as 
they  are  considered  as  a cause  of  the  significant  water  movements 
in  the  upper  layers  of  the  sea.  The  ocean  currents  are  generated 
and  maintained  chiefly  by  the  driving  forces  of  the  wind.  The 
"thermodynamic  engine"  (including  the^halinodynamic”)  is  relatively 
insignificant  in  the  oceans,  at  least  in  the  upper  layer?,  where 
currents  with  measurable  velocities  occur.  Weak  thermo-ha line  cir- 
culations, however,  will  dominate  in  the  deeper  strata  of  the 
oceans  beneath  a layer  which  will  be  called  the  "layer  of  no  motion" 
for  the  wind  induced  ocean  circulation  in  this  and  the  following 
reports. 

This  report  deals  only  with  the  wind  driven  ocean  circulation. 
The  temperature  and  salinity  structure  is  considered  only  as  far 
as  the  density  stratification  is  involved  in  the  final  development 
of  the  wind  driven  current  systems.  This  density  stratification 
is— beyond  the  minor  influence  of  "density  currents" — very  important, 
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since  in  a stratified  ocean  a ’’zero  layer"  (layer  of  no  motion) 
for  the  wind  driven  circulation  must  develop,  which  may  be  con- 
sidered as  the  lower  boundary  of  the  circulation,  where  the  velo- 
city of  the  wind  driven  currents  vanishes.  Therefore,  the  "zero 
layer"  can  also  be  replaced  by  a rigid  but  frictionless  boundary, 
and  the  necessary  resistance  in  order  to  balance  the  driving 
forces  is  brought  about  only  by  the  virtual  internal  friction. 

The  wind  driven  circulation  in  an  ocean  with  a variable 
<iar>t:h.  T).  n-f*  thfl  Mirrants  is  analysed,  and  it  is  found  that  the 

r 1 7 • — » • 

relative  change  cf  the  depth  D with  latitude  equals  the  relative 
change  of  the  Coriolis  parameter  in  the  large  scale  planetary 
circulation.  This  suggests  that  the  ocean  reacts  to  the  planetary 
vorticity  effect  in  such  a way  that  it  rather  tends  to  adjust  its 
level  of  no  motion  for  the  wind  driven  circulation  than  to  dis- 
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the  west.  Since  the  depth  of  the  layer  of  no  motion  for  the  wind 
driven  circulation  is  also  strongly  related  to  the  vertical  den- 
sity stratification,  the  mechanism  of  mutual  adjustments  between 
the  current  systems  ana  the  mass  distribution,  as  it  is  finally 
observed  in  the  oceans,  appears  rather  complicated. 

For  a quantitative  study  of  the  wind  driven  ocean  circulation 
Knowledge  of  the  wind  stress  at  the  sea  surface,  and  of  the  virtual 
Internal  friction  is  necessary.  Both  of  these  important  forces 
have  been  discussed  in  detail  in  this  report.  With  regard  to  the 
internal  friction  a simple  expression  corresponding  to  Guldberg- 
Mohn’s  expression  was  used,  and  coefficients  of  virtual  internal 
friction  wore  determined.  It  seems  that  a more  general  treatment 
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of  th9  internal  friction  will  not  essentially  change  the  results. 

The  opinions  of  different  authors  are  most  controversial  in 
discussion  of  the  wind  stress  at  the  sea  surface,  especially  at 
wind  velocities,  less  than  10  m/sec.  It  is  this  critical  range  of 
wind  speeds  which  comes  into  consideration  when  dealing  with  clima- 
tological wind  data,  over  most  parts  of  the  sea  surface.  At  present, 
more  and  more  evidence  has  been  accumulated  which  shows  that  the 
vljl yci iv6  wind  stress  st  low  wind  speeds  is  larger  then  it 
hitherto  assumed  by  aost  of  the  workers  on  this  subject.  There- 
fore, the  stress  rornruia  as  suggested  by  the  author  in  19*3  was 
used  to  calculate  the  wind  stress  for  given  wind  velocities. 

Some  special  cases  of  the  wind  driven  circulation  beneath  the 
layer  of  frictional  influence  are  analyzed,  such  as  the  meridional 
flow  along  the  east  coast  of  the  United  States,  and  the  zonal  flow 
of  the  Antarctic  Circum  Polar  Current.  The  outline  of  future 
work  briefly  discusses  the  ideas  which  may  lead  to  a better  approach 
of  wind  driven  circulation  in  a real  ocean,  taking  the  true  average 
wind  field  and  the  boundaries  into  account.  This  approach  is  based 
on  a numerical  integration  of  ths  differential  equation  for  the 
horizontal  mass  transport  in  an  ocean  where  the  wind  induced  cur- 
rent systems  have  a variable  depth. 
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NOTES  ON  THE  WIND-DRIVEN  OCEAN  CIRCULATION 


I.  Introduction 

Ocean  currents  with  the  greatest  velocity  are  observed  in 
the  upper  layer  of  the  sea  with  depths  varying  between  two  hundred 
and  two  thousand  meters.  This  circulation  is  generated  and  main- 
tained mainly  by  the  driving  force  of  the  wind.  The  question  is, 
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of  the  hydrodynamical  theories  as  they  have  been  developed  up  to 
the  present  time? 

The  problem  may  be  divided  into  two  parts.  Part  one  is  con- 
cerned with  the  "gross"  or  general  features  of  the  oceanic*  circu- 
lation, and  part  two  with  the  details  which  recently  have  been  re- 
vealed in  some  individual  branches  of  the  general  circulation  by  means 
of  refined  techniques  for  oceanographical  observations. 

There  is  no  doubt  that  not  only  the  details  but  also  the  gross 
features  of  the  oceanic  circulation  are  the  result  of  a complicated 
mechanism  in  which  different  forces  or  factors  are  involved.  It  is 
possible  that  the  mechanism  which  governs  some  outstanding  "gross" 
features,  such  as  the  "intensification"  of  the  Gulf  Stream  along 
the  east  coast  of  North  America  and  the  "intensification"  of  the 
Kuro  Shio  along  the  coast  of  eastern  Asia,  is  essentially  the  same 
as  the  mechanism  which  may  account  for  some  "details"  as  observed 
in  the  western  North  Atlantic  and  in  other  Individual  branches  of 
the  general  current  systems  of  the  oceans.  Such  details  are,  for 
example,  the  series  of  overlapping  currents  "arranged  somewhat 
like  the  shingles  on  a roof'  (Iselin,  1952),  and  certain  "counter 
currents."  It  seems  to  be  the  tendency  of  ocean  currents  in  strati- 
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fled  water  to  split  up  into  narrow  streams  of  higher  velocifcy 
rather  than  to  form  a broader  and  weaker  band  of  flow.  This 
"streakiness”  cf  ocean  currents  is  not  only  present  within  the 
swift  Gulf  Stream  and  Kuro  Shio,  but  also  in  other  branches  of 
the  circulation  along  the  eastern  and  western  boundaries  of  the 
oceans  as  well  as  in  mid=ocean  regions,  although  the  individual 
"streaks”  are  much  weaker.  At  present,  no  satisfactory  explanation 
for  such  details  in  the  flow  pattern  of  ocean  currents  has  been 


given. 

The  general  circulations  in  the  Atlantic  and  Pacific  Ocean 
are  asymmetrical  with  respect  to  the  equator,  and  this  asymmetry 
is  related  to  the  asymmetry  of  the  general  atmospheric  circulation. 
As  a consequence  of  this  a special  equatorial  current  system  is 
observed,  which  from  the  dynamical  point  of  view  need3  some  further 
consideration.  Common  to  the  oceans,  including  the  Indian  Ocean, 
is  the  fact  that  the  equatorial  currents  «re  relatively  strong 
compared  with  the  currents  of  the  subpolar  gyre.  This  fact,  as 
well  as  the  east-west  asymmetries  in  the  North  Atlantic,  in  the 
North  Pacific  and  partly  in  the  Indian  Ocean  cannot  be  explained 
in  a simple  way  by  the  wind  field  itself.  For  example,  the  Gulf 
Stream  and  Brazil  Current  are  almost  equivalent  branches  of  the 
subtropical  gyre  on  the  Northern  and  Southern  Hemisphere,  yet  the 
Gulf  Stream  transports  about  ten  times  as  much  as  the  Brazil  Cur- 
rent, ana  the  gross  features  ox  tue  equivalent  circulation  systems 
on  both  hemispheres  show  striking  differences. 

The  present  report  is  concerned  with  some  notes  on  the  general 
or  "gross"  features  of  the  wind  driven  circulation.  Although  no 
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final  answer  to  the  general  problem  of  the  horizontal  and  verti- 
cal current  stratification  can  be  given,  it  is  hoped  that  these 
notes  may  serve  to  stimulate  further  interest  in  the  subject 
along  lines  which  will  result  in  a resumption  of  the  problem  in 
the  way  it  was  initiated  by  Walfried  Ekman. 


II.  Historical  review 

Approximately  50  years  have  passed  since  Walfried  Ekman 
placed  our  conception  of  the  wind  driven  ocean  currents  on  a com- 
pletely new  basis.  In  his  well  known  theory  he  showed  how  the 
effect  of  the  wind  goes  far  beyond  the  generation  of  pure  drift 
currents,  and  how  the  gradient  currents  in  the  deeper  layers  are 
maintained  by  the  indirect  action  of  the  wind.  His  studies  were 
directed  toward  the  calculation  of  the  total  current  system  for 
the  given  wind  system  over  the  oceans.  He  tried  to  solve  the 
upoblss  f qt*  f ho  vsrt  icsl  snd  ths  hori vp!  o»a1  tv  d i s * ? i but  i or* « 


The  difficulties  which  rendered  Ekman 's  general  analysis  very  dif- 
ficult have  not  yet  been  overcome.  Ekman,  however,  succeeded  in 
solving  essential  partial  problems,  and  most  of  his  results  of 
the  first  paper  on  the  subject  (1905)  have  survived  half  a cen- 
tury; they  still  belong  to  our  basic  knowledge  in  dynamical  oceano- 
graphy. 

Ekman  was  the  first  to  introduce  eddy  viscosity  (virtual 
friction)  into  the  theory  of  ocean  currents.  In  a paper  from 
he  deduced  some  important  laws  which,  in  the  case  of  steady  motion, 
the  currents  and  particularly  their  curl  must  obey.  The  acceler- 
ation of  the  water,  originally  regarded  as  effectively  eliminated 
from  the  equations,  was  later  found  to  be  important  as  far  as  the 
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"deep-current " is  concerned.  In  his  revised  theory  of  1932*  the 
dynamical  theory  of  the  wind  driven  ocean  circulation  and  its 
problems  were  developed  by  Ekman  to  such  a degree  of  completeness* 
that  relatively  few  additions  have  been  made  since  that  time.  Prob- 
ably, the  only  point  where  Ekman* s model  of  steady  state  ocean 
currents  has  been  revised  essentially,  and  where  in  the  future 
essential  revisions  have  to  be  expected,  is  concerned  with  the 
mechanism  and  the  effect  of  internal  friction  and  vertical  density 
stratification. 

A new  epoch  of  theoretical  research  in  the  field  of  wind 
induced  ocean  currents  started  with  the  work  by  II. U.  Sverdrup 
(1947),  H.  Stommel  (1948),  K.  Fidaka  (1949)  and  W.  Munk  (1950). 

These  authors  essentially  gave  up  trying  to  evaluate  the  vertical 
velocity  distribution  and  confined  the  analysis  to  the  horizontal 
mass  transport  by  introducing  the  vertically  integrated  mass  trans- 
port as  the  dependent  variable.  This  was  necessary  in  order  to 
avoid  certain  difficulties  which  rendered  Ekman's  analysis  too 
complicated.  It  was  possible  to  examins  the  more  general  case  ox 
a baroclinic  ocean  without  having  to  specify  the  nature  of  the 
vertical  distributions  of  density  and  currents,  although  even  in 
this  case,  some  simplifying  assumptions  had  to  be  made.  The  prob- 
lem was  to  derive  the  gross  features  of  the  horizontal  mass  trans- 
port of  the  oceanic  circulation  for  a given  wind  system.  These 
stimulating  theories  can  be  considered  as  a remarkable  step  in 
dynamic  oceanography,  and  the  results  obtained  by  the  various  models 
are  of  great  interest. 

H.  Stommel  succeeded  in  1948  in  showing  that  due  to  the  effect 
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of  the  varying  Coriolis  force  with  latitude  the  gyre  of  the 
general  circulation  around  the  Sargasso  Sea  must  be  strongly 
displaced  toward  the  west  in  the  case  of  a constant  depth  of 
the  wind  driven  circulation.  The  results  are  very  strong  cur- 
rents along  the  east  coasts  of  the  continents.  This  '‘westward 
intensification"  of  ocean  currents  is  obvious  in  the  Gulf  Stream 
and  the  Kuro  Shio,  and  it  seems  that  a satisfactory  explanation 
of  this  striking  gross  feature  of  the  oceanic  circulation  has  been 
given.  However,  remarkable  exceptions  are  found  in  the  South  At- 
lont.in  ° nd  Scut*  Pacific , and  Stommel's  model  needs  fur-*^’*  con- 
sideration. 

Using  interesting  mathematical  methods,  Munk  and  collaborators, 
Hidaka  and  others  extended  the  theoretical  research  to  more  natural 
wind  distributions  and  more  realistic  ocean  boundaries. 

Hunk's  Investigation  (1950)  has  served  to  emphasize  the  fun- 
damental importance  of  the  wind  stress  curl  rather  than  of  the  wind 
stress  vector  in  determining  the  transport  of  ocean  currents  in 
meridionally  bounded  oceans.  In  the  case  of  an  irrotational  wind 
field,  such  as  a uniform  west  wlna  over  an  ocean  of  constant  depth 
bounded  on  its  western  and  eastgrn  sides,  it  follows  that  there  is 
no  net  mass t ransport , and  that  the  stress  is  balanced  by  the  pres- 
sure distribution  resulting  from  the  piling  up  of  water  against  the 
lee  shore.  This  result  becomes  understandable  only  when  keeping 
the  two-dimensional  model  and  its  theoretical  basis  in  mind.  The 
theoretical  basis  of  the  model  is  essentially  a differential  equa- 
tion representing  a vertically  integrated  version  of  the  vorticity 
equation  in  the  case  of  an  ocean  of  constant  depth.  It  expresses 
a balance  between  the  lateral  stress  torque,  the  planetary  vorticity* 
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and  the  wind  curl.  It  is  evident  that  the  vertical  velocity  com- 
ponent must  be  taken  into  account  in  order  to  explain  the  steady 
state  circulation  in  a bounded  ocean.  Otherwise,  the  result  that 
the  permanent  ocean  currents  are  only  related  to  the  rotational 

component  of  the  wind  stress  field  over  the  ocean  would  indicate 

/ 

that  the  currents  must  vanish  were  the  wind  stress  field  irro- 
tatlonal. 

Most  of  the  important  assumptions  of  the  previous  models  are 
concerned  with  the  lower  boundary  of  the  wind-induced  circulation 
system,  or  the  "depth"  of  the  currents,  and  with  the  concept  of 
lateral  and  vertical  eddy  viscosity,  especially  with  the  concept 
of  isotropic  lateral  mass  exchange. 

With  the  attempt  to  explain  some  of  the  observed  detail  flow 
patterns  in  the  Gulf  Stream  north-northwest  of  the  Azores,  in  th9 
California  current  and  the  Norwegian  current  by  taking  non-iso- 
t roplc  lateral  mass  exchange  into  account,  the  author's  attention 
was  called  to  some  probable  relationships  between  the  vertical 
density  gradient,  the  depth  of  the  current,  the  variation  of  lati- 
tude while  the  current  flows  in  meridional  direction,  and  the  "age" 
of  the  current,  that  is,  the  time  when  a wind  of  a certain  velocity 
starts  to  blow  over  an  ocean  originally  at  rest. 

The  varying  depth  of  the  wind  generated  oceanic  circulation 
is  the  first  fact  with  which  we  are  concerned,  and  which  has  been 
neglected  in  the  previous  models  even  in  the  case  of  a steady 
oceanic  circulation.  This  depth  depends  on  the  density  strati- 
fication of  the  ocean,  the  "age"  of  the  current,  and  possibly  on 
internal  friction  (lateral  and  vertical).  At  this  point  some 
difficulties  arose,  so  that  we  are  facing  the  problems  of  the  gross 
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features  of  the  steady  oceanic  circulation  again. 

In  the  present  report  it  will  be  shown  that  the  depth  of  the 
wind  induced  oceanic  circulation  plays  a significant  role  in  the 
hydrodynamic  analysis  of  the  problem.  A more  realistic  model  of 
the  wind  driven  circulation  of  the  North  Atlantic  Ocean  will  be 
treated  by  numerical  integration  of  the  basic  differential  equa- 
tions in  a subsequent  report. 

Ill*  The  differential  equation  of  horizontal  mass  transport 

Rectangular  coordinates  may  be  assumed  for  a first  simplified 
model.  Let  the  y-axis  point  northward,  the  x-axis  eastward,  and 
the  z-axis  vertically  upwards  from  the  undisturbed  sea  surface  which 
coincides  with  a level  surface.  The  boundaries  of  the  ocean  model 
may  be  considered  as  vertical  walls,  which  at  first  approximation 
would  account  for  the  continental  slope.  Besides  these  lateral 
boundaries  the  lower  boundary  of  the  current  system  has  to  be  fixed. 

It  is  known  that  in  equatorial  regions  the  layer  of  no  motion 
("zero  layer")  for  the  wind  driven  circulation  is  found  in  about 
200=300  m of  depth,  and  it  increases  with  increasing  latitude  to 
approximately  IJjOO  m in  45°-50°N.  (The  depth  of  this  "zero  layer" 
should  not  be  confused  with  the  depth  of  the  layer  of  frictional 
influence.)  There  is  much  evidence  for  this  fact,  and  it  has  been 
assumed  by  most  oceanographers  that  there  is  a continuous  change 
of  this  "level  of  no  motion"  between  low  and  high  latitudes.  Much 
support  of  this  point  of  view  has  been  given  by  A.  Defant  (1941) 
in  his  extensive  analysis  of  the  Atlantic  circulation,  where  inter- 
esting details  of  the  layer  of  no  motion  have  been  revealed. 

The  "zero  layer"  of  the  oceanic  circulation  in  the  wind  affected 
surface  and  subsurface  layers  may  be  considered  as  the  lower  boundary 
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of  the  current  system.  It  may  be  replaced  by  a rigid  surface 
where  the  velocity  vanishes  (no  friction  at  this  boundary).  There- 
fore, the  "zero  layer"  has  to  be  considered  as  the  depth  D of  our 
model  of  wind  driven  ocean  currents ? and  D is  a function  of  x and  y. 

Let  f = 2c»  sin9  be  the  Coriolis  parameter,  p the  pressure, 
u and  v the  horizontal  components  of  the  current  velocity  vector 
in  the  x and  y direction,  respectively,  v [cm”'1  gr  sec”1']  the  co- 
efficient of  eddy  viscosity  in  the  lateral,  and  fj.  the  coefficient 
of  eddy  viscosity  in  the  vertical  direction,  V the  Laplace  oper- 
ator in  the  x-y  plane  and  p the  density.  The  hydrodynamic  equa- 
tions of  steady  horizontal  motion  are 


pfv  + /x  + V V2u  - = 0 

► 

-pfu+/i.^— + vV2v  = o, 


one  hydrostatic  equation 
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and  the  equation  of  continuity 
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The  system  of  hydrodynamic  equations  (1)  to  (3)  together  with 
adequate  boundary  conditions  provides  the  basis  for  a very  general 
solution,  that  is,  to  find  the  horizontal  and  the  vertical  distri- 
bution of  u and  v.  Although  this  solution  is  of  great  interest, 
the  mathematical  analysis  is  too  complicated.  Therefore  the 
question  of  the  horizontal  mass  transport  may  be  considered  only. 
If  the  vertically  integrated  mass  transport  M,  with  the  components 
U,  V is  introduced  as  the  dependent  variable, 


U 


■L 


pu  da  , 


■/ 

-D 


ov  dz 


Since  the  integration  is  carried  out  from  a depth  where  the  total 
motion  is  zero  (u  = v = w = 0),  to  the  surface,  (z  - 0)*,  it  follows 
from  (3)  that 

. /°  r° 

f pu  dz  + 3^  / pv  dz  = 0 , 

-D  -D 

and  that 


dl™  * Si  * H = 0 


(4) 


Equation  (4)  states  that  the  horizontal 


iivergenoe  of  the 


total  mass  transport  equals  zero,,  in  contrast  to  the  divergence 
of  the  mass  transport  in  individual  levels.  For  example  in  the 
layer  of  frictional  influence  (wind  drift  current),  the  divergence 
may  be  zero  in  exceptional  cases  only. 


The  special  form  of  the  equation  of  continuity  for  the  verti- 
cally integrated  mass  transport  (4)  permits  the  introduction  of  a 
"stream  function",  v|/  , for  the  mass  transport,  and 

U = d \I//dy 

(5) 

v = - <>  \fj/dx 

Integration  of  equations  (1)  over  the  depth  between  the  level  of 


no  motion  (z  = -D)  and  the  surface  (z  = 0)  gives 


•"Actually  z = £ , where  £ the  elevation  or  depression  of  the  sea 
surface  with  reference  to  „the  level  surface  of  the  undisturbed  sea 
surface.  However,  since  Q is  a very  small  quantity  it  can  be 
neglected  against  D. 
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where  <}p/dx  and  dp/Oy  represent  average  values  of  the  horizontal 
pressure  gradients  in  the  whole  layer  between  z - -D  and  z = 0. 

The  two  integrals  in  equations  (6)  are  equal  to  the  differences 
of  the  components  of  the  tangential  stress  at  the  surface  and  at 
the  depth  -D.  The  tangential  stress  at  the  surface  can  be  evalu- 
ated from  the  wind  at  "anemometer  height"  by  means  of  an  empirical 
formula  (Neumann,  1948).  This  question,  however,  needs  some  further 
consideration  and  will  be  discussed  in  more  detail  in  section  IX 
of  this  report.  The  tangential  stress  in  the  depth  z = -D  equals 


wm  /\ 


Besides  the  wind  stress  at  the  surface,  horizontal  and  vertical 
shearing  stresses  in  the  water  have  to  be  taken  into  accounts  The 
dynamical  effect  of  these  internal  stresses  is  evident  in  the  vir- 
tual internal  friction.  Since  it  is  difficult  to  introduce  internal 
friction  in  a most  general  form,  and  in  order  to  simplify  the  dif- 
ferential aquations  for  the  mathematical  analysis,  it  seems  justi- 
fiable for  the  purpose  of  the  present  report  to  assume  that  the 
frictional  forces  are  proportional  to  the  total  mass  transport  M. 

The  combined  effect  of  the  virtual  friction  in  horizontal  and  verti- 


^ a -L 


A M « a4*  4 A 
Uli  CV/UJ, 


»m 

rxai  uc 


VI 4 0 u.  V 


~ ~ ~ ^ II  ~ 4 -~~!l 

a o au  ci,  x c v*  c#  x v c 


A ~4-  n 
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For  the  whole  water  column  between  z = -D  and  z = 0 the  components 
of  the  effective  internal  frictional  forces  are  assumed  to  be 


» , . r dL 

“x  dv 


Rv  = + r 


s*-.J 


(7) 


where  r[sec“^]  represents  an  effective  "friction"-coef f icient . A 
plausible  average  value  of  this  coefficient  has  to  be  determined 


later. 
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With  these  assumptions  equations  Co)  can  be  written 


£- 


T + r 
x 


T - r 

y 


sf -iD-»i 


T and  r are  the  components  of  the  wind  stress  at  the  sea  sur- 

^ y 

face  in  the  direction  of  the  x-axis  and  y-axis,  respectively* 

Cross  differentiation  of  (8),  ana  subtraction  leads  to 

V Jt  U.Y 

According  to  equations  (8) 

/)—  % rx  An  in 

v V ou  VJK  2^  _ 

<Jx  oy  “ 5y  ox 

By  substituting  this  into  equation  (9)-  it  follows 

- * w 

- a Tx  a Ty  Tx  dD  Tx  dD  ni^ 

~5y  " <3x  D oy  D ox 

With  this  differential  equation  and  adequate  boundary  conditions 
the  horizontal  mass  transport  is  uniquely  determined  in  the  ocean 
under  consideration.  The  boundary  conditions  for  the  North  Atlantic 

Qt^aan  all  1 ha  spec  if  i e d IstSf. 

In  the  special  case  where  the  depth  D is  assumed  to  be  constant, 


it  follows  from  (11) 


l n Itlf 


= - curlz  T , for  D = const. 

This  equation  is  identical  with  Munk's  differential  equation  for 
the  horizontal  mass  transport  (1950$  equ.  6),  if  the  sign  of  p is 
changed  (Munk  defined  V = d P/Cz,  U = - dP/dy)  and  the  term  rV2  p 
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is  replaced  by  in  Munk's  notation. 

Before  an  evaluation  of  equation  (11)  several  important  factors 
involved  in  the  model  have  to  be  discussed  and  to  be  fixed.  This 
is  concerned  with  D as  a function  of  x and  y,  with  the  wind  stress 
at  the  sea  surface  and  with  the  effective  friction  coefficient, 
r,  in  the  water  as  defined  before. 

IV.  The  " zero-layer"  in  the  Atlantic  Ocean 

A successful  attempt  to  determine  the  depth  of  the  zero  layer 
(layer  of  no  motion)  in  the  Atlantic  Ocean  was  made  by  A.  Defant 
(1941).  When  comparing  the  differences  in  the  dynamic  depth  of 
given  pressure  surfaces  between  adjacent  oceanographic  stations, 
it  was  found  that  in  certain  levels  this  difference  was  almost 
constant  over  a rather  large  depth  interval.  This  would  mean  that 
either  the  whole  body  of  water  throughout  this  depth  interval 
(sometimes  with  a thickness  of  500-800  meters)  has  a constant  velo- 
city in  the  vertical  direction,  or  that  thi3  layer  is  absolutely 
at  rest. 

Defant's  analysis  of  the  North  and  South  Atlantic  Ocean  has 
shown  that  it  is  much  more  plausible  to  assume  that  this  layer  is 
absolutely  at  rest,  than  that  it  is  in  relative  motion  to  some 

- A.  1 n ^ •» *»  Jl  4 - J.  1 T\ — ^ 4. J 4.V.- 

uoilttl*  ioioAdhwd  jlcvox  m ouo  uvujr  ui  Warn'#  wx  cui aoouuww,  uua 

average  depth  of  this  layer  to  be  the  "level  of  no  motion"  for  his 
dynamical  computations.  This  hypothetical  level  of  no  motion  was 
traced  throughout  the  whole  North  and  South  Atlantic,  and  it  was 
found  that  it  composes  without  any  constraint  to  a closed  layer 
of  varying  thickness.  This  fact  and  others  as  pointed  out  by 
Defant  may  serve  to  strengthen  our  confidence  in  Defant’s  analysis 
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of  the  circulation  of  the  Atlantic » Any  other  assumption  of  the 
level  of  no  motion  would  lead  to  unreasonable  and  contradictory 
results  about  the  vertical  structure  of  the  Atlantic  circulation. 

The  same  hypothesis  was  applied  to  a dynamical  analysis  of  the 
Black  Sea  circulation  by  the  author  (194-2,1943).  In  this  case 
the  results  could  be  checked  against  other  techniques  in  deter- 
mining the  level  of  no  motion,  and  the  computed  dynamical  struct  tire 
of  the  Black  Ssa  was  in  good  agreement  with  the  observed  conditions* 

Defant 's  results  for  the  Atlantic  Ocean  are  shown  in  figure  l. 
The  depth  of  the  "zero  layer"  (D)  increases  with  increasing  lati- 
tude, both,  on  the  Northern  and  on  the  Southern  Hemisphere.  In  the 
North  Atlantic,  the  depth  of  the  zero  layer  is  more  irregular  than 
in  the  South  Atlantic,  which  has  its  analogy  in  the  more  compli- 
cated current  pattern  of  the  North  Atlantic. 

Based  on  the  chart  shown  In  figure  1,  zonal  averages  of  the 
depth  of  no  motion  (D)  were  computed  for  the  North  Atlantic,  and 
plotted  against  the  latitude  in  figure  2.  For  the  South  Atlantic, 
individual  values  ware  used  as  given  along  the  meridian  of  20°W, 
since  the  variation  of  D in  the  east -west  direction  is  small  com- 
pared with  the  variation  of  D in  the  meridional  direction.  The 
values  of  D,  according  to  A.  Defant's  chart,  are  indicated  in  figure 
2 by  circles,  whereas  the  curve  represents  the  function  D = -K  sinq>. 
The  constant  K is  different  on  the  Northern  and  Southern  Hemisphere, 
but  the  variation  of  D with  latitude  follows  this  law  closely.  It 
is  seen  that  the  zonal  averages  of  the  depth  of  the  layer  of  no 
motion  in  the  North  Atlantic,  and  the  corresponding  depths  along 
20®W  in  the  South  Atlantic,  can  be  represented  very  closely  by 

D(cp)  = -K  sincp,  (12) 
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Fig*  ...  Depth  of  the  "Zero  Layer"  in  the  Atlantic  Ocean, 

according  to  A*  Defant  (19lil).  The  figures  at  the 
lines  are  to  be  multiplied  by  100,  thus  denoting 
the  depth  in  meters. 
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Figure  2.  Average  depth,  D,  of  the  layer  of  no  motion  in  the  Atlantic  Ocean. 
The  curves  « — represent  D (?)  ■ const.  sin<jp  , and  the 

values  according  to  A.  Defant  (I9I4I)  are  indicated  by  circles  (0). 
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or,  the  relative  variation  of  D is  given  by 

i = Jf  ctn  n>  , 

1/  yy  « 

where  R is  the  radius  of  the  earth.  Since 

t = R otD  f > 

there  is  much  evidence  that  in  the  "gross"  features  of  the  oceanic 
circulation 

dT  _ f 3D 

3y  = 5 5y  ’ VXJ' 

and  the  important  term 

\ V#Jf  i-S  v/  ' VA 

in  equation  (11)  drops  out  in  the  "gross”  features  of  a model  where 
only  zonal  conditions  are  considered. 

In  the  model  by  Stommel  (1948),  Hidaka  (1949)  and  Munk  (1950) 
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D was  assumed  to  be  constant.  Thererore,  the  term  responsible 
for  the  "westward  intensification",  , which  contains  the 

variation  of  the  Coriolis  parameter  with  latitude  (planetary  vor- 
ticlty)  was  not  cancelled  out  by  the  variation  of  the  depth  D with 
latitude  (topographic  vorticity,  according  to  Ekraan's  terminology). 

It  seems  unlikely  that  the  east-west  asymmetry  of  the  gross 
features  of  the  oceanic  circulation  is  caused  by  the  variation  of 
the  Cor  lulls  parameter  with  latitude,  but  since  there  Jjs  an  inten- 
sification of  ocean  currents  along  the  east  coasts  of  North  America 
and  Eastern  Asia,  this  particular  flow  pattern  requires  another 
explanation.  This  may  also  answer  the  question  why  an  analogous 
intensification  of  the  subtropical  gyre  (Munk,  1950)  is  not  found 
in  the  Southern  Hemisphere. 


It  seems  more  likely  that  the  ocean  reacts  to  the  planetary 
vorticity  effect  in  such  a way  that  it  rather  tends  to  adjust  its 
level  of  no  motion  for  the  wind  driven  circulation  than  to  displace 
the  whole  gyre  with  relatively  high  current  velocities  to  the  west . 
Whether  or  not  this  adjustment  is  complete  is  a question  wnich  still 
has  to  be  answered. 


It  may  be  added  that  the  relationship  as  stated  in  equation  (13) 
does  not  say  anything  about  "cause  and  effect".  The  zero  layer  nay 
be  the  result  of  the  current  system,  and  probably  it  is  to  some 
degree  in  the  gross  features  and  even  in  the  details.  However, 
several  other  external  and  internal  forces  and  factors,  besides  those 
taken  into  account,  are  involved  in  the  mutual  relationships  between 
the  oceanic  circulation  and  the  stratification.  Some  of  these 
factors  are  the  large  scale  heat  gain  and  loss  through  the  sea  sur- 
face and  evaporation,  precipitation,  and  lateral  and  vertical  mixing 
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processes.  The  first  mentioned  climatic  factors  and  mixing  processes 
are  especially  responsible  for  some  large  scale  density  strati- 
fications. With  "large  scale  density  stratification"  in  this  sense 
is  meant,  for  example,  the  very  weak  vertical  density  gradient  in 
polar  regions  and  the  strong,  significant  vertical  density  gradient 
in  tropical  regions,  both  taken  as  an  average  between  the  surface 
and  the  bottom  water.  The  vertical  density  stratification  is  an 
important  factor  in  the  oceanic  circulation,  which,  at  least  more 
indirectly,  may  determine  the  wind  driven  circulation.  This  will 
be  shown  by  a simple  model  in  the  next  two  sections. 


V . A simple  model  of  a stratified  ocean 

Disregard  at  first  the  differences  in  the  salinity  distribution 
and  assume  the  temperature  distribution  as  the  main  factor  of  the 
density  stratification.  The  vertical  temperature  distribution  as 
observed  in  the  deeper  layers  of  the  oceans  is  to  a certain  degree 

~ A s*  « « 4-  U ^ /vb  4-1 — J ~ -4 ..1  ~ 4-  _ 1 

uux  xjr  ou  a q uiic  xuouiv  oi  vna  uv^oa^ia.^  uxi  vuaaux^ii  ui  one? 

"spreading"  of  water  masses  (G.  Wust,  1935)*  This  is  only  partly 
true,  as  shown  by  H.U.  Sverdrup  (1939).  Stationary  conditions  which 
appear  on  the  average  "gross"  features  to  be  rather  similar  to  the 
observed  conditions  in  the  ocean  may  exist  even  in  the  case  where 
the  components  u,  v,  w,  of  the  oceanic  circulation  are  zero;  that 
is,  in  an  ocean  completely  at  rest. 

Assume  in  the  first  step  an  ocean  at  rest.  This  "ocean"  may 
cover  the  whole  earth,  and  only  the  sun's  radiation  (with  zero  de- 
clination) may  warm  up  the  equatorial  regions.  The  depth  of  the 
ocean  may  be  constant  (=  h meters).  In  this  case,  the  resulting 
temperature  distribution  will  be  symmetrical  in  the  meridional  di- 
rection and  equal  in  the  zonal  direction.  Since 


u - v = w = 0 by 


assumption, 


8 - fe  £ 8)  ♦ 8>  ♦ Ms1 


P ox / o y 'p  oy'  os  'P  oz 
where  T is  the  temperature,  and  Ax , Ay,  Az  are  exchange  coefficients 
for  heat  conduction.  With  an  east-west  symmetry,  as  in  our  model, 
we  need  only  to  consider  the  resulting  temperature  distribution 
in  a meridional  cross  section.  Stationary  conditions  are  possible 


if 


^ fy)  = - h (^*$1)  • 


Due  to  the  sun's  radiation  there  is  a heat  gain  at  the  equator 
and  a heat  loss  at  the  poles.  Suppose  that  a meridional  temperature 
distribution  is  maintained  at  the  sea  surface,  and  let  T(y)  at  the 
sea  surface  be 

T(y)  = Tp  + To  003  §5  y , (14) 

where  Tp  is  the  temperature  of  the  sea  surface  at  the  pole, 

(y  = b),  Tp  + Tq  the  temperature  at  the  equator  (y  = 0). 

If,  for  simplicity,  we  assume  *y  and  A„  to  be  constant  in 
space,  but  Ay>-  Az,  a solution  of  the  differential  equation 

(15) 


Ay  + Az  = 0 
* q y o z 


with  the  boundary  condition  (14)  is 

T(y,z)  = Tp  + T0  ^°^o°hhah  ^ ’ 005  lb  y • 

where 


(16) 


a = 

' *D  ✓ A- 


/5 


Assume  that  the  ratio  Ay/Az  - 10',  and  then,  with  b = 90°  - 10,000  Km, 

a = 0.497  • 

Further,  let  T0  = 2?eC,  Tm  = 0°C,  and  let  h = 6000  m. 


IS 
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With  the  assumptions  made,  stationary  conditions  of  the 


temperature  distribution  in  an  ocean  at  rest  ?culd  ho  pvooxulo  xi 


the  temperature  in  different  depths  at  different  latitudes  follows 
equation  (16).  In  the  following  table  this  temperature  distribution 
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is  computed  on  the  assumption  that  Ay/A^  = 10' 


This  temperature  stratification  (figure  3a)  shows  rough  simi- 
larities with  the  observed  conditions  in  its  very  gross  features. 
These  similarities  are  the  strong  vertical  stratification  in  tropical 
and  subtropical  regions  with  decreasing  vertical  gradients  towards 
the  bottom  and  the  homogeneity  at  the  poles.  The  distribution  of 
the  isotherms  in  this  model  (figure  3a)  may  lead  to  the  impression 
that  the  cold  water  from  polar  regions  spreads  out  along  the  ocean 
bottom  towards  equatorial  regions.  Even  a more  realistic  picture 
could  easily  be  derived  with  the  assumption  of  an  asymmetrical  tem- 
perature distribution  at  the  Surface  (Northern  Hemisphere  warmer 
than  the  Southern). 

VI.  Notes  on  the  Influence  of  the  wind 

In  the  second  step,  suppose  that  the  wind  blows  over  this  ocean. 
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Fig.  3a.  Comp” ted  temperature  distribution  in  an  ocean  at  rest  (equilibrium 
of  heat  conduction  with  given  sea  surface  conditions  Ay/Az  - 10f), 

3b«  Adjustment  of  mass  distribution  (temperature)  to  wind  induced  cii> 
culation  with  the  additional  assumption  of  geos trophic  equilibrium 
(schematic).  The  line  -+-  -4-  -4-  indicates  the  layer  of  no  motion 
for  the  wind  induced  circulation. 
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The  prevailing  winds  are  east  winds  between  the  equator  and  30°, 
west  winds  between  30°  an,d  60°,  and  weaker  east  winds  between  60° 
and  the  pole.  The  primary  effect  of  the  winds  acting  at  the  sea 
surface  is  to  produce  drift  currents  with  velocities  which  decrease 
with  depth  and  which  are  essentially  limited  to  the  upper  50  m or 
100  m of  depth  (layer  of  frictional  influence).  In  our  model,  fig- 
ure  3b,  which  is  considered  to  represent  the  Northern  Hemisphere, 
the  total  mass  transport  of  water  in  the  layer  of  frictional  in- 
fluence (d)  is  directed  toward  the  north  between  the  equator  and 
30°N,  and  between  60°N  and  the  pole,  whereas  the  total  mass  trans- 
port is  directed  towards  the  south  between  60°  and  30°.  The  con- 
sequences are  a divergent  flow  in  equatorial  regions  and  around  60° 
and  a convergent  flow  in  the  latitudes  around  30°  and  near  the  pole. 
The  divergence  (or  convergence)  of  this  mass  transport  is  propor- 
tional to  the  curl  of  the  wind  stress,  and  it  can  easily  be  computed. 

The  unequal  divergence  at  different  localities,  which  strongly 
depends  on  details  of  the  wind  field,  sets  up  a slope  of  the  sea 
surface,  in  our  simplified  model,  which  drops  from  30°  toward  the 
equator  and  towards  60°  and  rises  again  from  60°  toward  the  pole  as 
indicated  schematically  in  figure  3*>.  From  this,  hydrostatic  pres- 
sure gradients  result  in  the  water,  and,  as  an  indirect  consequence 
of  the  wind  effect,  gradient  currents  will  develop  in  the  deeper 
layers  (in  addition  to  the  "density  gradient  currents",  simply  called 
"density  currents",  due  to  the  original  horizontal  temperature  dif- 
ferences). The  depth  of  the  wind  induced  gradient  currents  depends 
on  the  vertical  stratification  of  the  ocean  water. 

In  what  follows,  we  may  neglect  the  "density  currents"  as  they 
would  result  from  the  original  temperature  distribution  in  the  model. 
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To  examine  qualitatively  the  wind  effect  on  the  stratified  ocean, 
as  shown  in  figure  3a,  let  us  assume  that  the  velocity  of  the  wind 
induced  gradient  currents  decreases  with  increasing  depth.  In  the 
steady  stat9  the  sea  surface  slopes  must  have  reached  a certain 
stationary  inclination,  and  the  mass  distribution  (temperature  field) 
must  have  changed  in  such  a way  that  beneath  the  layer  of  frictional 
influence  from  point  to  point  the  geostrophic  equilibrium  condition 


is  nearly  fulfilled,  or  that,  approximately, 

du  _ g 1 dp 

<7z  “ f p oy 

holds.  Here,  u is  the  zonal  velocity  of  the  gradient  current  in 
the  x-direction  (indicated  in  figure  3b  by  feathered  arrows),  p the 
density  (which  in  our  model  is  only  a function  of  the  temperature 
at  constant  pressure),  and  y the  inclination  angle  between  the  lines 
of  equal  density  and  the  (horizontal)  level  surfaces.  Qualitatively, 
the  distribution  of  isotherms  in  the  meridional  cross-section  under 
these  conditions  should  look  like  the  dashed  lines,  as  shown  in 
figure  3*>* 

When  looking  in  the  direction  of  the  current  between  60°N  and 
30°N  the  isotherms  must  slope  down  from  the  left  hand  side  to  the 
right  hand  side,  that  is,  from  north  to  south,  provided  the  current 
decreases  with  increasing  depth,  and  they  must  slope  down  from  the 
equator  towards  30°N,  since  the  wind  induced  gradient  current  in 
this  region  flows  from  east  to  west.  This  leads  to  a thick  warm 
water  accumulation  around  30oN9  and  to  a relative  cooling  in  corres- 
ponding depths  in  equatorial  regions.  Furthermore,  the  vertical 
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temperature  gradient  in  certain  subsurface  levels  between  30°  and 
the  equator  must  iflcreass,  and  beneath  this  region  decrease , when 
compared  with  the  original  conditions  for  an  ocean  at  rest,  since 
the  geostrophic  equilibrium  condition  requires  a tilting  up  of  the 
isotherms  toward  the  equator  in  the  layer  of  the  equatorial  (grad- 
ient) current.  The  depth  of  this  layer  depends  on  the  stratification, 
and  it  is  indicated  by  the  line  d — k-4-  in  figure  3b. 

VII«  The  hydrostatic  mass  compensation 

The  condition  of  a steady  state  requires  that  for  each  column 
of  water,  the  divergence  of  the  total  mass  transport  M from  the  stir- 
face  to  the  bottom  or  to  the  level  of  no  motion  equals  zero.  Let 
M-^  be  the  mass  transport  of  the  drift  current  (in  the  layer  of 
frictional  influence),  and  let  be  the  mass  transport  of  the 
gradient  current  (from  the  surface  to  a depth  where  the  current 
vanishes);  then  div  (M-,  + 1^)  = 0,  or  div  M-,  = - div  1L,  for  each 

WU  «—  A ■ 

column.  Only  the  total  mass  transport,  integrated  from  the  surface 
to  the  layer  of  no  motion,  is  non-divergent . But  the  layer  of 
frictional  influence  is  a divergent  layer,  and  this  divergence  must 
be  balanced  by  divergent  gradient  currents  in  deeper  layers. 

In  the  case  of  an  ocean  with  constant,  or  almost  constant- 
potential  density  from  the  surface  to  the  bottom,  the  hydrostatic 
pressure  gradient,  as  caused  by  the  "wind  set  up",  should  practical- 
ly extend  down  to  the  bottom  (as  in  figure  3b  near  90°  of  latitude). 
Here,  a "bottom  current"  as  assumed  in  Ekman’s  "elementary  current 
system"  would  provide  the  necessary  anisobaric  mass  transport  in 
order  to  balance  the  anisobaric  mass  transport  ir«  the  layer  of 
frictional  influence.  However,  in  a strongly  stratified  ocean,  (as 
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in  fig.  3b  between  30°  and  0°),  a layer  of  no  motion  may  develop 
(with  steady  state  conditions)  in  a certain  depth  above  the  bottom 
by  •'mass-compensation”  of  the  hydrostatic  pressure  differences  re- 
sulting from  sea  surface  slopes.  In  this  case,  the  necessary  balance 
for  the  mass  transport  in  the  layer  of  frictional  influence  has  to 
be  provided  for  completely  by  the  wind  induced  gradient  current. 

If  the  hydrostatic  equation  holds  in  vertical  direction,  the 
pressure  p at  a depth,  h,  is  given  by 

z=0 

P = - g / P <iz  + g p0£0 
z=-h 

where  is  the  elevation  of  the  see  surface  above  the  level  sur- 
face z = 0,  and  pQ  is  the  density  cf  the  surface. 

Assume  two  oceanographic  stations,  A and  B,  fairly  close  to- 
gether, for  which  the  vertical  density  stratification  is  known,  and 
suppose  L £ is  the  elevation  of  the  sea  surface  at  station  B re- 
lative to  station  A where  CA  =0.  Then,  the  depth  of  the  level 
of  no  motion  for  the  current  perpendicular  to  the  line  AB,  between 
these  stations  is  given  at  an  average  (constant)  depth,  z = -D,  where 


and 


u 


. r = 


. JL 

Po 


f 


(o . 


A 


-D 


Or,)  dz  . 

* n 


(17) 


Let  the  coordinate  system  be  oriented  with  the  x-axis  in  the  (zonal) 
wind  direction.  Then,  in  the  layer  of  frictional  influence 
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(18) 


(Ml)y  = " 

according  to  Ekraan. 


x 

f * 


‘Vx  = 0 


The  wind  induced  gradient  current  is  a non-divergent  current 

only  as  a first  approximation.  As  an  effect  of  internal  friction 

there  must  be  a slight  deviation  from  the  direction  parallel  to 

the  isobars.  The  average  angle  of  deviation  (from  high  to  low 

pressure)  in  the  gradient  layer  may  be  denoted  by  o.  Let  u be 

6 

the  zonal  component,  of  the  geostrophic  current  in  the  case  of  no 


internal  friction,  then 


= . 1 tE 

f dv  » 


Dg  = / dz  * - f / dz 

-D  -D 

Suppose  that  the  actual  total  horizontal  mass  transport,  U,  in 
the  gradient  current  is  deflected  slightly  from  the  direction  of 
U by  the  angle  a,  as  schematically  indicated  in  figure  4.  The 

b 

magnitude  of  a was  estimated  as  l°-2°  in  higher  latitudes  and 
5°-6®  in  lower  latitudes  (Neumann,  1952)* 


Low 


Hiah 


**)  — ♦ 

Fig.  t.  The  difference  vector  U - represents  approximately  the  aniso- 
baric  mass  transport  to  balance  the  mass  transport  in  the  layer 
of  frictional  influence. 


The  difference  vector  ff  - U represents  approximately  the 

o 

aniscbaric  mass  transport  Mg  to  balance  the  mass  transport 
in  the  steady  state,  and 


It,  = - i sine 


(19) 


With  the  hydrostatic  equation 

0 


and 


= ' 8 / 
-D 


dp 

ft 


dz 


/ 5y  dz  * - e/fZ5y  d2j 

-D  -D  J 


dz 


the  mass  transport  Mg  can  be  written,  approximately, 


-*  ■ h 


dz 


dz 


(20) 


for  two  closely  spaced  stations,  A and  B,  with  a distance  A y apart. 

From  the  steady  state  condition,  div  M = 0,  and  from  ( 1 8)  and 
(20)  we  have  with  the  assumption  that  between  the  two  closely 
spaced  stations  the  variation  of  and  Mg  can  be  neglected 


T 


X 


(21) 


This  equation  can  be  solved  numerically  for  D,  the  level  of  no 
motion.  The  zonal  wind  stress,  Tx,  can  be  determined  from  the 
wind  field,  and  with  regard  to  c the  assumption  may  be  made  that  a 
in  the  region  of  the  equatorial  current  is  approximately  5°-6°  and 
decreases  to  i.5c  in  60s  latitude. 
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Figure  5 shows  the  zonal  average  density  distribution  in  a 
meridional  cross  section  through  the  South  Atlantic  Ocean.  The 
zonal  average  of  the  density  was  computed  from  the  data  as  given 
by  G.  Wust  and  A.  Defant  (1936)  in  the  atlas  of  the  "Meteor”  work# 
From  Schott’s  (194-2)  wind  charts,  zonal  average  wind  velo- 
cities in  the  South  Atlantic  have  teen  derived,  and  the  wind  stress, 
Tx,  has  been  computed  by  means  of  an  empirical  formula  (Neumann, 
1948).  More  information  about  the  wind  stress  is  given  in  section 
IX  of  this  report.  The  zonal  wind  velocities  W„,  the  wind  stress 
T^,  and  the  estimated  values  for  a are  given  in  Table  1. 

Let  us  assume  equally  spaced  "stations"  with  a meridional 
distance  of  A y - 5°j  between  5«5°S  and  57«5°S.  From  (21)  with 
Ay  = 555  km  and  <Tt  = (p  - 1)  x 10^,  we  have  for  an  estimate 

dz  (22) 


The  average  vertical  density  distribution  as  shewn  in  figure 
5 is  tabulated  for  the  numerical  integration  in  table  2. 

Table  1 


Latitude 

A 

V 

t- 

10 

15 

20 

25 

Wx(m/sec) 

-0.5 

-3.0 

>•4.0 

-5.3 

-3.5 

0 

Tx( dyne/cm2) 

0.05 

0.57 

0.9 

1.38 

0.74 

0 

a° 

- 

__7.7 

6.2 

4.7 

3.9 

3.2 

Latitude 
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40 

45 
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Q.A 
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•7.4 
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Tx(dyne/cnr) 

0.50 

1.22 

2.2  7 

3.25 

2.87 

2.36 

0.97 

a0 

2.7 

2.3 

2.1 

1.9 

1.7 

1.6 

1.5 
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Latitude 


Average  meridional  density  section  (Of)  through  the  South  Atlantic  Ocean,  and  average  depth 
of  the  layer  of  no  motion  according  to  A.  Defant  (— t—  — t-  -*-)  compared  with  calculations  (®) 


Table 


For  two  fictitious  "stations"  at  42.5°3 
ively,  the  following  vertical  distribution  of 


and  4-7. 5°S,  respect 


I 


OV  >dz 


dz 


is  derived 

• 

depth 

Aert  * 

(m)  * a 

- % 

• Az 

Z(*V  Az) 

T 

TAz 

I 

x 104 

x 104 

x 104 

x 108 

x 108 

0 

0.40 

0 

0 

0.64 

0.32 

0.64 

0.64 

200 

0.24 

0.64 

A Mr* 

0.84 

1.68 

400 

0.16 

1.04 

2.32 

600 

0.33 

1.205 

2.410 

0.17 

1.37 

4.73 

800 

0.34 

1.54 

3.08 

0.17 

1.71 

7.81 

0.16 

0.33 

1.875 

3 .75 

1000 

2.04 

11.56 

1250 

0.375 

2.228 

5.57 

0.14 

2.415 

17.13 

A AA 

V/#  J \J\J 

2.565 

6.42 

1500 

0.10 

2.715 

23.55 

0.08 

0.225 

2.828 

7.08 

1750 

2.940 

30.63 

0.137 

3.085 

7.70 

2000 

0.03 

3.077 

38.33 

2500 

0.100 

3.127 

15.64 

0.01 

3.177 

53.97 

0.025 

3.189 

15.95 

3000 

0.00 

3.202 

69.92 

With 

= 3.25 

A 

dyne/cm  , and  a = 1.9°, 

Tx  . 
sina 

5.65  X 

107  = 

53  x 108, 


therefore  D = ?50Q  a for  the  zonal  component.  This  compares  fairly 
well  with  the  depth  of  the  layer  of  no  motion  as  shown  in  figure  2. 
Similar  computations  for  higher  latitudes  (55°S)  lead  to  values 
which  are  a little  too  large,  whereas  the  depth  D computed  for  the 
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South  East  Trade  wind  region  (10°S  through  35°S)  agree  well  with 
Defant 's  values.  The  computed  values  D are  platted  in  the  den- 
sity section,  figure  for  comparison  with  the  depth  D in  the 
central  parts  of  the  South  Atlantic  according  to  figure  2. 

Probably,  the  estimated  angles  a are  the  most  inaccurate. 
Hcwever,  the  estimate  of  D by  equation  (22)  should  serve  only  to 
check  the  order  of  magnitude  of  the  facto*  j involved  in  this  pos- 
sible mechanism  of  •’mass  compensation”  to  determine  the  level  of 
no  motion.  Another  fact  which  must  be  kept  in  mind  is  that  this 
calculation  only  considers  the  meridional  component  of  the  ani- 
sobarlc  mass  transport. 

At  this  point,  the  dashed  curve  in  figure  5 may  at  least  be 
mentioned  briefly.  It  connects  the  maximum  depth  of  the  lines  of 
equal  density  ( 0*t) , and  is  in  some  way  related  to  the  tilting  of 
the  axis  of  the  subtropical  gyre.  This  tilting  in  meridional 
direction  toward  the  poles  with  increasing  depth  has  already  been 
mentioned  by  A.  Defant  (1941),  but  so  far  no  explanation  has  been 
given  for  this  gross  feature  of  the  oceanic  circulation. 
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The  coefficient  of  internal  friction 

In  figure  6,  M is  the  vector  of  the 
© 


J II 

geostrophic  mass  trans- 


port, and  M is  the  vector  of  the  true  mass  transport  which  is  de- 
flected by  the  angle  a from  the  direction  parallel  to  the  isobars. 

The  difference  vector,  - M , is  in  some  way  related  to  the  average 

w 

internal  friction  in  the  whole  water  column  from  z = 0 to  z = -D. 

Assume  horizontal,  non-accelerated  motion  under  the  influence 
of  internal  friction.  Then 

1 ^ -4  -4  *4 

- “ grad  p+M  xf+RyO 
P 


> 


(23) 
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Figure  6.  Force  diagram  of  a non-accelerated  gradient 
current  with  friction. 


whore  R is  the  vector  of  the  internal  frictional  force.  For  a 


geostrophic  current  without  friction. 

- “ grad  p + Mg  x f = 0 
On  substituting  (24)  into  (23),  we  have 

(M  4)  x f = R . 

O 


(24) 


(25) 


Thus,  the  vector  of  internal  friction  is  perpendicular  to  the 

->  — > 

difference  vector  M - M,  and  only  in  exceptional  cases  is  p equal 

o 

to  zero. 

Suppose  a is  known,  then  p and  I R I = R can  be  computed.  In 
equation  (8),  for  the  components  of  the  effective  internal  friction, 
we  used  for  a first  approximation  Guldfeerg=Mohn'3  expression 


Rx  = - rU 
Ry  « - rV  , 


by  assuming  R 5 - rM.  As  in  figure  6,  let  G be  the  vector  of  the 


i 

) 


average  pressure  gradient  force  and  C the  vector  of  the  Coriolis 


force. 


and 


O XUV^v 


Then  from  figure  6, 

G sino  = rM  cosp  (26) 

G cosa  ” 2o>  sin<p  M + rM  sinp  , (2?) 

the  components  of  the  forces  in  the  direction  of  M and  per- 


pendicular to  it  must  balance.  M - M is  proportional  to  the 

o 

average  '‘dragging  force",  according  to  Sprung. 

From  (26^  and  (27) 


ct  ns 


2 m sincp 
r cosp 


♦ tanp 


(28) 


2®  sin®  M 
tanp  = etna  - -g-^- 

In  the  case  where  p -*  0 

- 2m  sin® 

etna  = „ — *•  . 

r 

So  far,  no  direct  computations  of  r have  been  made  in  dynamic 
oceanography,  and  it  seems  that  the  only  estimates  of  a are  those 
used  in  this  report-  (Neumann.  1 952 ) . 

Whether  (30)  can  be  used  instead  of  (28)  to  calculate  r from 
a,  depends  on  the  order  of  magnitude  of  p.  However,  data  which  are 
accurate  enough  to  compute  p are  not  available. 

The  mass  transport  of  the  Antarctic  Clrcum  Polar  Current  between 
65°S  and  45°3,  is  given  by 

M » 1014  [gr  sec"1]  , 

approximately,  according  to  Sverdrup  et  al  (194-6).  (This  value  may 
be  even  a little  larger.)  The  average  meridional  pressure  gradient 
G can  be  estimated  fairly  well  from  the  "Meteor"  data  published  by 
G.  Wust  (1939).  Using  "Meteor"  Stations  Mt  130,  124,  60,  and  10, 
the  following  pressure  values,  p,  in  decibar  are  obtained,  relative 
to  the  sea  surface  where  p = 0? 


(29) 

(30) 
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p_  in  decibars 

- z 

Depth  (m) 

Mt.  130 
q>  = 63?7ftS 
Ar,  laiS- 

Kt.  324 

<P  - 55- s 
X - 16 

Mt.  60 
® = 48.5°S 
A = 16,9°W 

Mt.  10 
a = 41,2°S 

XsLiZiA3 

0 

0 

0 

0 

0 

50 

51.368 

51.361 

51.358 

51.336 

100 

102.772 

102.740 

102.730 

102.682 

150 

154.192 

154.137 

154.114 

154.043 

200 

205.626 

205.552 

205.513 

205.425 

300 

308.529 

308.428 

308.355 

308.237 

400 

411.481 

411.361 

411.252 

411.110 

6oo 

617.534 

617.388 

617.209 

617.021 

800 

323.783 

823.616 

823.377 

823.137 

1000 

1030.228 

1030.039 

1029.758 

1029.461 

1200 

1236. 869 

1236,660 

1236.343 

1236.000 

1400 

1443.703 

1443.47 7 

1443.127 

1442.745 

l6CO 

1650.728 

1650.489 

1650.109 

1649.695 

2000 

2065.352 

2065.092 

2064.659 

(2064.191) 

2500 

2584.70 

2584.41 

2583.93 

- 

3000 

3105.24 

3104.92 

3104.39 

- 

From 

these  data  G 

- 1,4  X 10" 

4 J ,3 

dyne/cnH . 

With  M « 1014  gr 

sec"^*,  D s 

= 3000  m and  1 

a total  width  of  2.22  x 

103  Km,  it  follows 

per  unit  1 

area  of  the  cross  section 

that 

M[gr 

-1  -2.-1 
sec  cm  J - 

1014  gr 

sec"3- 

s 1.5  , 

2.22  x 3.0  : 

TO  O ~ 

x 10--’  cm 

and,  with  a - 1 .7°,  according  to  Table  1,  for  tp  = 55% 

tanp  s 34  - 43  . 


This  result  can  only  serve  to  check  the  order  of  magnitude  of 
a,  G,  and  M,  since  tanp  is  given  as  a small  difference  between 
two  larger  quantities  which  are  known  only  approximately.  There- 
fore, p can  not  be  determined  from  these  data.  However,  since  the 
deviation  of  the  true  mass  transport  from  the  direction  of  the 
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geostrophic  mass  transport  is  very  small;  the  assumption  may  be 
justified  that  p ~ 0 for  the  purpose  Of  estimating  r.  Then,  from 
(30)  and  table  1 for  <p  = 55* 

r - = 3.34  x IQ*6  [sec"1]  . 


etna  - — -1 

In  the  lower  layers  of  the  atmosphere,  r is  on  the  average 
approximately  100  times  larger.  H.U.  Sverdrup  and  J.  Holtsmark 
(1917)  derived  for  10  North  American  stations  average  values  from 
50-100  individual  observations.  They  found  r between  7 x 10  ' and 
2.7  x 10“4  sec"'1',  and  p between  -12®  and  +48°. 

The  height  of  the  anemometer  above  the  ground  is  of  great  in= 
fluence.  Both  r and  p decrease  with  height,  as  shown  by  Sverdrup 
and  Holtsmark: 


height  (meters) 

25 

25-49 

50 

v x 104 

2.20 

1.99 

1.56  sec 

P 

zy 

20- 

rr 

With  increasing  height,  p and  r decrease  in  the  lower  layer  of 
the  atmosphere,  but  it  seems  that  no  fairly  well  established  in- 
formation about  these  friction  coefficients  is  available  for  the 
higher  troposphere,  which  is  comparable  with  the  layer  of  wind 
driven  ocean  circulation  beneath  the  layer  of  frictional  influence. 
It  may  well  be,  that  in  the  higher  troposphere  p— *-0,  magn.  a s 2°; 
and  magn.  r « 10  ' - 10  . 

According  to  equation  (30)  etna  should  vary  approximately 
proportionally  to  sir. 9 , if  r is  considered  to  be  constant.  Since 
a is  given  in  table  1 by  previous  results,  r in  different  lati- 
tudes is  shown  in  the  following  table: 
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Latitude 

10° 

0 

0 

01 

_30° 

40° 

50« 

60  c 

U 

6,2 

3.9 

2.7 

2.1 

1.7 

1.5 

r x 1C6 

2.75 

3.14 

3.43 

3.44 

3.30 

3.30 

With  the 

estimated  angles 

a,  r is 

almost 

constant.  This  may 

indicate  that  in  the  average  for  the  whole  layer  of  the  gradient 
current,  frictional  conditions  are  rather  uniform  in  the  different 
branches  of  the  circulation  system,  although  under  special  con- 
ditions significant  local  variations  of  the  effective  frictional 
forces  may  occur.  For  the  numerical  integration  of  the  differential 
equation  (11),  however,  an  average  value 

r « 3,3  i 10“^  [sec“'L] 

will  be  used. 

IX.  The  wind  stress  at  the  sea  surface 

Estimates  of  the  total  wind  stress  at  the  sea  surface  rest 
Inrgelv  on  observation??  of*  the  vert  ion!  distribution  of  mean  wind 
near  the  sea  surface.  Such  wind  profile  measurements  In  the  lowest 
layers  over  the  sea  surface  have  been  compared  with  the  velocity- 
distribution  in  circular  pipes,  so  that  the  interfacial  stress  could 
be  deduced  from  the  measured  wall  stress  in  pipes.  This  method 
simply  applies  the  results  of  experimental,  aerodynamical  studies, 
and  is,  therefore,  based  on  the  assumption  that  such  a comparison 
is  justified. 

Regardless  of  the  question  as  to  how  far  the  analogy  between 
the  wall  stress  in  pipes  and  the  wind  stress  at  the  free  sea  surface 
goes,  there  are  still  other  objections  against  the  method  used  , 
especially  with  regard  to  the  method  of  observing  the  average  verti- 
cal wind  profile  with  cup  anemometers  at  fixed  ’’heights55  above  the 
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undisturbed  sea  level  (or  above  the  highest  wave  crests). 

Most  of  the  observations  seem  to  indicate  that  with  an  adia- 
batic lapse  rate  in  the  lower  layers  above  the  sea  surface  the 
observed  wind  profile  can  be  approximated  fairly  well  by  a loga- 
rithmic law  within  a certain  height  interval.  However,  in  the  lowest 
layers,  some  50  or  100  cm  above  the  crests  of  the  sea  surface  rough- 
ness elements  (waves),  the  logarithmic  curve  gradually  steepens 
more  and  more  when  approaching  the  "sea  surface”.  This  feature  is 
characteristic  for  the  observed  wind  profiles  over  the  rough  sea 
surface,  and  a careful  analysis  of  the  observations  shows  that  the 
logarithmic  law  for  the  wind  profile  Just  above  the  waves  is  a 
rather  questionable  approximation. 

The  value  and  accuracy  of  wind  profile  measurements  with  cup 
anemometers  for  evaluating  the  wind  stress  at  the  sea  surface  have 
been  discussed  in  some  detail  by  the  author  (1949,  195D-  There 
are  mainly  two  facts  which  make  this  method  unfit  for  stress  com- 
putations: 

First  of  all,  the  questions  remains  open  as  to  what  level  has 
to  be  chosen  as  the  "height  of  the  sea  surface",  or,  in  other  words, 
what  is  the  "sea  surface"  in  the  case  of  an  air-sea  interface  with 
traveling  waves?  Any  attempt  to  derive  a vertical  wind  profile 
from  observed  wind  speeds  at  fixed  heights  above  the  crests,  re- 
quires a clarification  of  this  question. 

Secondly,  wind  measurements  with  cup  anemometers  at  fixed 
heights  above  the  wave  crests  are  unfit,  to  derive  true  average  wind 
speeds  over  the  rough  sea  surface.  The  wind  at  a fixed  height  over 
the  trough  of  a "wave"  may  have  only  40%  - 5 0 % or  less  of  its  speed 
over  the  crest.  With  steep  waves,  even  a "lee -eddy"  may  develop 
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somewhere  near  the  trough,  and  the  wind  direction  may  be  reversed. 
Since  wind  profile  measurements  near  the  surface  have  been  taken 
only  at  wind  speeds  of,  approximately,  less  than  10  m/sec,  the 
rather  short  period  wave  motion  passing  by  underneath  the  cup  ane- 
mometers "speeds  up"  the  anemometers  every  time  a crest;  passes. 

The  Inertia  of  the  cnemometers  does  not  &?„lcw  slowing  down  over 
the  troughs,  and  therefore,  almost  maximum  wind  speeds  are  observed 
Instead  of  average  values.  When  used  for  stress  computations,  too 
small  values  of  the  stress  are  obviously  derived  by  this  procedure. 
The  "critical  wind  speed"  at  7 m/sec  as  discussed  by  W;  Munk  (1947) 
can  be  explained  as  a result  of  \.he  inaccuracy  of  the  data  obtained 
by  wind  profile  measurements  (Neumann,  195D  • 

Another  more  direct  method  of  deriving  the  wind  stress  at  the 
sea  surface  is  based  on  observations  of  the  steady  state  sea  sur- 
face slope  in  which  the  water  is  piled  up  due  to  wind  ("wind  set 
up").  This  method  is  the  oldest,  and  it  was  first  used  by  W.  Ekman 
in  1905*  Since  that  time  much  more  data  on  the  'Wind  sat  up"  in 
different  seas  and  lakes  have  accumulated,  and  with  these  data,  the 
method  was  reexamined  and  critically  applied  by  the  author  In  1948. 

The  stress,  T , may  be  expressed  in  terms  of  a resistance  co- 
2 

efficient,  y , which  is  defined  by 

T = P*  Y2  W2  , (31) 

where  p’  is  the  density  of  the  air,  and  W is  the  mean  wind  speed  at 
a given  distance  above  the  sea  surface.  In  practice*  this  distance 
is  given  by  the  "anemometer  height",  about  10  m.  The  question  is 
whether  y is  a constant,  or  whether  it  depends  on  the  wind  speed 
itself.  In  the  case  of  a flow  over  a solid  surface  with  well  de- 
fined and  constant  roughness  conditions,  such  as  grass  land  or  snow, 
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Y*  may  be  a constant  for  different  wind  speeds.  However,  in  the 

case  of  the  sea  surface,  the  "roughness"  i3  given  by  the  composite 

wave  motion.  Since  the  wave  motion  is  a function  of  the  wind  speed, 

W,  (and  of  the  state  of  wave  development),  y will  probably  also  be 

a function  of  the  wind  speed.  (In  the  case  of  stable  stratified 

air,  y also  depends  on  the  stability  of  the  air  stratification.) 

The  results  obtained  by  the  author  in  1948  indicate  that  under 

normal  conditions  (with  an  adiabatic  lapse  rate  in  the  lowest  layers) 
2 

Y varies  inversely  proportional  to  the  square  root  of  W,  and 


..2  _ „ „ _ 3 / 1 [cm/sec]  \ 1/2 

t 7.^  a j.  \ yi  j^cm/sec]^ 


(32) 


for  winds  from  1 to  40  m/sec-  Thus*  the  wind  stress  at  the  sea 


surface  (for  a fully  arisen  sea)  can  be  represented  by 


T = p ’ x 


0.009  x W3/2  [dyne/cm2]  , 


(33) 


where  W is  to  be  inserted  in  centimeters  per  second. 


a cuiaparxsun  ux  \$<z)  wztn  tne  c orresponamg  resistance  co- 
efficients obtained  from  wind  profile  measurements  was  made  in 

the  results  are  reproduced  in  figure  y,  where  more  recent 
cup  anemometer  measurements  by  U.  Roll  (1948)  show  again  the  char- 
aeteristic  discrepancy  between  the  y values  as  derived  from  "slope" 
observations  and  from  wind  profile  observations.  By  disregarding 
the  results  based  on  "slope"  observations  at  wind  speeds  of  less 

than  10  m/sec,  and  replacing  them  by  the  results  of  the  inadequate 

2 

cup  anemometer  measurements,  a "Jump"  of  v from  low  to  high  values 
npnoar>  erroneously  at  10  m/sec  wind  velocity. 

A critical  review  of  the  accuracy  and  a summary  of  the  results 
obtained  by  different  methods,  including  wind  slope  observations 
in  wind  tunnels  has  been  recently  made  by  R.  B.  Montgomery  (19^2). 
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Fig.  7.  ttesistance  coefficiar.t  (y*  ~ 


■5")  at  the  sea  surface 


as  a function  of  wind  velocity  (W  m/sec),  according  to 
slope-observations  (wind  set  up)  and  wind  profile 
measurements*  (After  G.  Neumann,  1951). 


He  finally  comes  to  the  rather  disappointing  conclusion  that  at 
no  wind  speed  is  the  stress  (or  the  resistance  coefficient)  known 
confidently  to  within  half  its  value. 

At  present,  a new,  direct  approach  for  measuring  the  inter- 
facial stress  at  the  sea  surface  is  being  undertaken  by  L.A.E.  Doe 


as  part  of  this  project.  It  is  hoped  to  derive  the  surface  stress 
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from  direct  measurements  of  the  eddy  stress,  -p[w'u’]>  at  a short 
distance  from  the  interface,  and  from  wind  profile  measurements 
at  constant  heights  above  the  true  wavy  sea  surface  with  a special 
type  of  practically  inertialess  anemometers. 

For  a quantitative  study  of  the  wind  driven  oceanic  circulation, 
the  exact  knowledge  of  T (or  'f‘~)  is  very  important.  Since  the 
average  climatological  wind  data  used  in  such  a study  are  less  than 
10  m/ sec  in  most  areas  over  the  oceans,  the  most  contraversial  re- 
gion with  W < 10  riv/sec  comes  into  consideration  in  a critical  way. 

W.  iiunk  (1950)  in  his  study  of  the  wind  driven  oceanic  circulation 
came  to  the  conclusion  that  the  stress  values  as  derived  from  wind 
profile  measurements  are  much  too  small  to  supply  the  necessary 
driving  force  of  the  wind  at  low  wind  speeds.  The  author  agrees 
completely  with  W.  Munk's  statement  that  the  transports  of  the 
large  oceanic  currents  are  probably  as  good  an  indicator  of  the 

overall  stress  exerted  by  the  winds  on  the  ocean  as  any  of  the 

2 

measurements  on  which  the  value  of  y is  now  based.  Similar  con- 
clusions have  been  drawn  with  regard  to  the  wave  generating  forces 
by  wind  (Neumann,  1950). 

Recently,  in  a careful  analysis  of  storm  tide  data  on  Lake 
Erie,  G.  H.  Keulegsn  (1952)  also  derived  wind  stresses  for  larger 
wind  speeds  which  are  in  substantial  agreement  with  the  values  sug- 
gested by  the  author  (194-8). 

There  is,  at  present,  no  evidence  which  really  shows  dis- 
agreement with  the  use  of  the  empirical  stress  formula  as  proposed 
by  equation  (31)  with  the  variable  resistance  coefficient  of  (32). 
Therefore,  after  a critical  reexamination  of  the  papers  on  this 
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subject  since  1948,  this  formula  has  been  applied  for  computing 
the  wind  stress  in  table  1 of  this  report,  and  it  is  also  being 
used  for  a numerical  analysis  of  wind  driven  circulation  in  the 
North  Atlantic  Ocean. 


X.  Conclusions,  and  some  remarks  on  the  model  of  the  wind  driven 
ocean  circulation 

a)  The  differential  equation  (11)  for  the  horizontal  mass 
transport  of  the  wind  driven  oceanic  circulation  seems  to  promise 
a more  realistic  approach  to  a quantitative  analysis  of  the  hori- 
zontal mass  transport  in  the  upper  layers  of  the  sea.  It  comes 
a little  closer  to  Ekman's  original  formulation  of  the  problem, 
since  the  “topographic  curl  effect'1  was  shown  to  be  of  significant 
importance  throughout  all  of  his  papers. 

By  writing  equation  (9)  in  terms  of  the  vertically  Integrated 
mass  transport,  M,  with  the  components,  U and  V,  it  follows  that 

V A-f  1 / Ari  AT)  AyT  AT)\  -i 

curiz  M = eurV  , U3> 

where  dp/dx  and  dp/dy  again  represent  average  values  of  the  hori- 
zontal pressure  gradients  in  the  whole  layer  between  z = -D  and 


z = 0.  With  Ap 
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On  substituting  (34)  into  (33)  we  have 


curl. 

z 


and  with  M = 


- curl  R 


(35) 


Only,  if  D « constant 

- — d f , 

+ curiz  « = - v f curxz»  , v^o; 

and  the  equation  expresses  a balance  between  the  curl  of  internal 
stresses  (friction),  and  the  planetary  vorticity  and  the  curl  of 
the  wind  stress.  The  differential  equation  was  used  in  this  form 
in  the  preceding  papers  on  the  wind  driven  circulation. 


• v « a 1 1.1.  . ^ j -t-i.  1 — j .i  J n / n c'N  4 . 
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preted  as  the  z-component  of  a vector  product  of  two  vectors, 
grad  D = a and  grad  c -t  . If  6 is  the  angle  between  the 
vectors  a and  b , then 


as  _ aL  $2  = 

qx  gy  gy  gx 


ab  sind  = 0 


would  mean  that  a \\  b , and  the  lines  of  equal  dynamic  height  of 
the  sea  surface  (relative  to  a level  surface)  are  parallel  to  the 
lines  of  equal  depth  D* 

Similarly,  the  term 

HE  s 0 

dx  dy  “ <)y  dx 

would  express  the  requirement  that  the  lines  of  equal  D and  the 
isopycnals  p are  parallel.  Approximately  35  years  ago,  Bj„  Helland- 
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Hansen  and  W.  Ekman  (1923)  noticed  that  over  large  areas  in  the 
Norwegian  Sea  the  large  scale  spatial  distribution  of  temperature 
and  salinity  (density),  followed  closely  the  contour  lines  of  the 
bottom  topography.  Furthermore,  when  studying  large  scale  conditions 
in  the  oceans,  they  found  a remarkable  parallel  between  the  iso- 
therms and  isohalines  and  the  contour  lines  of  the  isobaric  surfaces . 
This  important  law  became  known  as  the  "law  of  the  parallel  sole- 
noids.” Although  this  law  seems  very  well  verified  in  the  gross 
features  of  the  large  scale  conditions,  there  are  slight  deviations 
which  are  nevertheless  very  important. 

Precise  leveling  along  the  American  East  Coast  shows  that  the 
mean  sea  level  increases  towards  the  north  from  Florida  to  Nova 
Scotia  (G.  Dietrich,  1937).  The  magnitude  of  this  slope  is  10“^, 
and  thus  It  is  of  the  same  order  as  the  sea  surface  slopes  derived 
from  data  in  the  Caribbean  Sea  (H.U.  Sverdrup,  et  al. , 1946).  In 
the  Caribbean  Sea,  as  in  many  other  parts  of  the  ocaa^,  the  slope 
of  the  sea  surface  is  practically  compensated  for  by  the  distri- 
bution of  mass.  However,  this  is  not  the  case  in  the  region  of 
the  Gulf  Stream  along  the  east  coast  of  North  America.  G*  Dietrich 
(1937)  has  shown  that  along  the  continental  slope  the  distribution 
of  density  docs  not  indicate  any  rise  of  the  sea  surface,  and  this 
is  true  also  in  the  case  where  the  sea  surface  tonoeranhv  is  re- 
f erred  to  the  .2000  dbar -surf ace  instead  of  che  oxygen  minimum  layer, 
which  was  used  by  Dietrich  (H.U.  Sverdrup  et  ai.,  1946). 

With  the  assumption  that  within  the  layer  of  the  wind  induced 
gradient  current  curl2  R S curl^'  , then  from  (35) 

v $-£-!§)  = <>. 


or 


v#-«[(Dll  + f||-)$f-(D3|+5||-)|i]  =0  (37) 

The  x,y-coordinate  system  can  be  oriented  in  such  a way  that 
the  y axis  points  in  a direction  i?  of  the  flow  of  the  Gulf  Stream 
along  the  east  coast  of  the  United  States.  Then,  the  £ -axis, 
perpendicular  to  ?}  , points  toward  the  southeast,  approximately. 

In  this  case,  it  is  found  from  figure  1,  that  du/d£  z -2.5  x 10“^ 
between  Cape  Hatteras  and  Nova  Scotia.  With  the  assumption  that 
dD/d7?  s 0,  or  at  least  dD/d7?  < dD/d£  , it  follows  from  (37)  that 

v $£  c°s  45°  + g(D  + p j^-)  = 0 . (38) 

v -f  v * i v •[  w- 

This  equation  relates  the  (average)  density  distribution  and  the 
slope  of  the  sea  surface  in  the  direction  of  the  flow  with  the  slope 
of  the  layer  of  no  motion  perpendicular  to  the  flow  and  the  verti- 
cally integrated  mass  transport  of  the  Gulf  Stream  while  moving 


northeastward. 


Suppose  that  there  is  no  compensation  of  di>  /d7?  by  the  density 
distribution  p {rj  ),  and,  therefore,  in  the  direction  of  flow  dp/d7? 

» 0.  Then  d£  /d7?  must  be  positive,  since  V(df/dy)  is  a positive 
quantity,  and  dD/d£  < 0,  The  order  of  magnitude  of  d£  /d7?  can 
be  estimated  with  the  assumptions  made  before,  and 


4 » iJLiZLl?!  . v | ctn? 

i gp  • 2 . 5 • 10  gp  • 2 . 5 • 10  J 


(39) 


between  Cape  Hatteras  and  Nova  Scotia  where  9 = 39°,  and  R = 6370  km 
is  the  radius  of  the  earth.  According  to  H.U.  Sverdrup  et  al.(I946) 
the  total  mass  transport  of  the  Gulf  Stream  towards  the  NE  equals 
55  Mill,  per  sec.  With  a total  width  of  the  current  of,  approxi- 
mately, 250  km, 


4=r 
1 ✓ 


and  from  (39) 


V » 2.2  x 10^  [gr  cm"1  sec**1], 


aC  /d7?  - 1.15  X 10“7  . 

This  compares  well  with  the  observed  magnitude  of  the  slope 
between  Cape  Hatteras  and  Nova  Scotia  (Sverdrup  et  al.,  1946)*  If 
this  slope  did  not  exist,  and  the  assumptions  on  which  (38)  is  based 
hold,  p would  have  to  increase  toward  the  NE  in  the  axi"  of  the  flow 
If  the  term  V(df/dy)  would  be  disregarded,  it  is  seen  from 
(38)  that  an  upward  slope  in  the  direction  of  flow  must  be  compen- 
sated for  by  a decrease  of  the  average  density  of  the  Gulf  Stream 
water  off  the  continental  shelf. 

It  should  be  kept  in  mind  that  equation  (38)  holds  only  in  the 
case  of  vanishing  by  small  friction.  It  certainly  does  not  apply 
to  the  surface  layers  of  frictional  influence,  and,  if  present, 
to  a bottom  layer  of  frictional  influence  (Ekman's  bottom  current). 
Similarly,  the  law  of  the  parallel  solenoids  does  nut  apply  to 
layers  where  frictional  influences  dominate. 

c)  With  a discussion  of  equation  (11)  for  the  horizontal  ma3s 
transport,  the  case  of  a simple  zonal  model  may  be  considered, 
where  (see  section  IV) 


_ I dD 
5y  D oy  > 


y v 0 y 

and  any  effects  of  meridional  boundaries  are  disregarded.  Further, 
suppose  that  T = T (y),  T = 0,  dD/dx  = 0,  V = 0.  With  these 

A A J 

assumptions , equation  (11)  reduces  to 

ilt  . 1 42  di.  = 1 ( 11s  . IiPi 

^y2  D dy  07  * Koy  D 5y<  ’ 


or 
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(40) 


By  integration,  it  follows  that 

rj  _ 45k.  = JLi  + c . D 
dy  r ’ 

where  C is  a constant,  which  becomes  zero  if  U = 0 for  T__  = 0* 

This  result  shows  that  the  zonal  mass  transport,  TJ,  is  proportional 
to  the  zonal  component  of  the  wind  sires*,  T^.  This  is  different 
from  Munk's  general  result  where  in  a confined  ocean  of  constant 
depth  it  was  found  that  the  horizontal  mass  transport  of  permanent 
ocean  currants  depends  only  on  the  rotational  component  of  the  wind 
stress  field  over  the  ocean.  Of  course,  in  an  ocean  bounded  by  con- 
tinents, tne  assumptions  v/hieh  led  to  (40)  are  in  general  not  justi- 
fied any  longer.-  Even  in  the  case  where  T may  be  assumed  to  be 

J 

zero,  &D/dx  is  not  necessarily  zero,  and  V ^ 0,  at  least  for  some 
distance  along  the  boundaries.  The  whole  system  of  horizontal  mass 


transport  will  be  modified  by  tne  boundary  conditions,  dug  De  sides 
a dependence  of  the  mass  transport  on  the  curl  of  T , there  will 
also  be  a dependence  on  T . 

Equation  (40),  however,  will  hold  approximately,  in  the  central 
parts  of  the  oceans,  far  away  from  the  boundaries.  Probably,  it 
can  best  be  applied  to  the  Antarctic  Circum  Polar  Current,  just  where 
W.  Munk  and  E.  Palmen  (1951)  found  discrepancies  with  the  results 
of  Munk's  theory. 


In  the  Antarctic  Circum  Folar  Current  between  65*3  and  45°S, 


according  to  table  1,  f » 2.0  dyne/cm^,  ana  r 
(see  page  30).  Thus  according  to  (40) 

U = 6 x 10^  gr  cnf^  sec”'*', 


3.34  x 10"“ 


sec 
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and  the  average  total  mass  transport  through  a section  between 
65°  and  45°S  is,  therefore, 

6 x 10*  x 2.22  x 10S  = 1,33  x 1014  [gr  sec-1], 

T 

which  agrees  well  with  the  observed  mass  transport.  According  to 
AoJ.  Clowes'  (1933)  computations  the  mass  transport  through  the 
Drake  Passage  is  about  1.1  x 10AM-  gr  sec  , and  it  is  very  probable 
that  the  transport  in  the  central  parts  of  the  South  Atlantic  Ocean 
is  approximately  20#  greater. 

XI.  Outline  of  future  work 

Equation  (11)  is  being  used  for  determining  the  horizontal 
mass  transport  of  the  wind  driven  circulation  in  the  North  Atlantic 
Ocean  between  the  equator  and  about  55°N.  The  east  and  west  bound- 
aries of  the  ocean  are  the  continents.  They  are  approximated  as 
closely  as  possible  by  broken  straight  lines  following  the  conti- 
nental shelf  along  the  200  m isobath. 

The  boundary  condition  along  the  continental  slopes,  which  may 
be  identified  with  the  continents  themselves,  requires  that  the 
mass  transport  be  parallel  to  the  boundary  curves.  Since  an  additive 
constant  with  the  stream  function  y makes  no  difference,  the  boundary 
condition  is  given  by  U/  = 0 along  the  continents. 

However,  since  this  model  is  not  completely  enclosed  by  land, 
more  complicated  conditions  have  to  be  introduced  along  the  northern 
and  southern  boundaries,  and  between  the  Antilles  Arc  and  Florida. 

The  Guiana  Current  along  the  northeast  coast  of  South  America  trans- 
ports an  appreciable  amount  of  water  across  the  equator  which  enters 
the  model  ocean  in  the  southwest  corner.  This  fact  cannot  be  dis- 
regarded. It  has  to  be  taken  into  account  as  a "source"  by  the 
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boundary  conditions  along  the  equator.  Another  "source",  which 
is  probably  very  important,  is  the  flow  of  wator  through  the  Straits 
of  Florida. 

Along  the  northern  boundary  a "source"  has  to  be  considered 
northeast  of  the  Grand  Banks  due  to  the  inflow  of  the  Labrador  Cur- 
rent. An  important  "sink"  between  Great  Britain  and  Iceland  is 
given  by  the  transport  of  the  North  Atlantic  Current  which  accounts 
for  a net  balance  of  the  water  in  the  region  covered  by  the  model. 

The  total  net  intake  of  water  between  the  equator  and  the  northern 
boundary,  that  is,  the  sum  of  the  "sources"  and  "sinks",  must  be 
zero. 

For  a numerical  integration  of  the  differential  equation  (11) 
the  surface  of  the  ocean  is  thought  to  be  covered  by  a series  of 
lattice  points  as  given  by  the  intersections  of  a set  of  rectangular 
lines  in  the  north-south  and  the  east-west  direction.  The  lattice 
of  this  net  consists  of  "interior"  points  and  "boundary-points". 
Interior  points  are  called  such  points  which  are  surrounded  by  ad- 
jacent points  belonging  to  the  lattice  of  the  covered  ocean  area. 

An  interior  point,  therefore,  is  always  surrounded  by  four  adjacent 
points,  whereas  a "boundary  point"  has  less  than  four  adjacent  points. 
The  unit  mesh  width  of  the  lattice  was  chosen  as  d = 222  km*  how- 
ever, in  the  western  North  Atlantic  it  may  be  necessary  to  use  a 
smaller  unit  (55*5  km)  in  order  to  get  some  of  the  details  of  the 
Gulf  Stream  flow  pattern.  This  section  will  be  considered  by  a 
special  evaluation  of  the  basic  differential  equation. 

The  differential  equation  (11)  for  the  horizontal  mass  trans- 
port can  be  written  in  the  form 
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rV2'!'  + FCz>y)  + G(x,y)  - H(x,y)  = 0 , 

whore  the  coefficients  F,  G and  H have  been  plotted  as  functions  of 
x and  y on  charts.  The  wind  stress  was  computed  for  the  month  of 
February  by  means  of  formula  (33) > using  the  wind  charts  put  out 
by  the  Marine  Branch  of  the  Meteorological  Office,  Air  Ministry, 
London  (19*8). 

If  the  subscripts  j,k  indicate  an  interior  point,  and  the  ad- 
jacent points  of  the  lattice  are  denoted  by  subscripts  (j,k-l), 
(j,k+l),  (j-l,k),  (j+l,k),  respectively,  the  difference  equation  for 
the  horizontal  mass  transport  becomes 

/ \lr  1 lr-1  + 1 1/1  W.1  + V + ]r  “ 4\fj  i 

r ^ r - r * I_>L:±1S L.lrA 1* L-ili  j + 

-3fk  * oj>k  . %k  = 0 . 

Such  a difference  equation  can  be  established  for  each  interior 
point.  There  are  as  many  equations  as  interior  points,  and  since 
the  number  of  interior  points  equals  the  number  of  unknowns,  with 
the  boundary  conditions,  \f/  Is  uniquely  determined  in  the  interior 
of  the  area,  provided  that  the  determinant  of  the  system  of  equations 
Is  unequal  to  zero. 


Basic  data  for  the  determination  of  the  coefficients  F,  G and 
H as  functions  of  x and  y?  and  tb*  results  of  the  numerical  .inte- 
gration will  be  published  in  the  next  report. 
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